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We study the effect of electron interactions in topological crystalline insulators (TCIs) protected 
by mirror symmetry, which are realized in the SnTe material class and host multi-valley Dirac 
fermion surface states. We find that interactions reduce the integer classification of noninteracting 
TCIs in three dimensions, indexed by the mirror Chern number, to a finite group Zg. In particular, 
we explicitly construct a microscopic interaction Hamiltonian to gap 8 flavors of Dirac fermions on 
the TCI surface, while preserving the mirror symmetry. Our construction builds on interacting edge 
states of Uil) x Z 2 symmetry-protected topological (SPT) phases of fermions in two dimensions, 
which we classify. Our work reveals a deep connection between 3D topological phases protected by 
spatial symmetries and 2D topological phases protected by internal symmetries. 
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The prediction and observation of topological crys¬ 
talline insulators (TCIs) in the SnTe material class has 
expanded the scope of topological matter and gained 
wide interest W- These TCIs possess topological sur¬ 
face states that are protected by mirror symmetry of the 
rocksalt crystal and become gapped under symmetry¬ 
breaking structural distortions [6H9]. These surface 
states are predicted to exhibit a plethora of novel phe¬ 
nomena ranging from large quantum anomalous Hall con¬ 
ductance [Biiniiii] to strain-induced pseudo-Landau lev¬ 
els and superconductivity [T2] , which are currently under 
intensive study [TSHIS]. 

According to band theory, TCIs protected by mirror 
symmetry are classified by an integer topological invari¬ 
ant, the mirror Chern number |16| . However, recent the¬ 
oretical breakthroughs [HHlll have found that the clas¬ 
sifications of interacting systems are markedly different 
from noninteracting systems in various classes of topolog¬ 
ical insulators and superconductors protected by inter¬ 
nal symmetries [25j . This raises the open question about 
the classification of interacting TCIs protected by spatial 
symmetries. On the experimental side, a growing body of 
interaction-driven phenomena has been found in existing 
TCI materials, including spontaneous surface structural 
transition and gap generation [SHS] and anomalous bulk 
band inversion [55]. Moreover, new TCI materials have 
been predicted in transition metal oxides HZl |5H] and 
heavy fermion compounds [291130) , where strong electron 
interactions are expected. 

Motivated by these theoretical and experimental de¬ 
velopments, in this work we study the effect of electron 
interactions in mirror-symmetric TCIs. Our main result 
is that interactions reduce the classihcation of 3D TCIs 
from Z in the noninteracting case to Zg. We obtain this 
result by introducing a “domain wall” construction of in¬ 
teracting surface states of 3D TCIs, which exploits the 
nonlocal nature of mirror symmetry. This construction 
builds on interacting edge states of 2D TCIs or 17(1) x Z 2 
symmetry-protected topological (SPT) phases, which we 


classify. Our work reveals a deep connection between 3D 
topological phases protected by spatial symmetries and 
2D topological phases protected by internal symmetries. 

Interacting TCIs in two dimensions: We first 
study interacting TCIs in two dimensions to set up the 
basis of later analysis in three dimensions. These 2D sys¬ 
tems have two independent symmetries: the 17(1) charge 
conservation and the mirror symmetry under the reflec¬ 
tion z —>■ —z, where z is normal to the 2D plane. Since 
this mirror symmetry is a Z 2 internal symmetry |31j . 2D 
TCIs with mirror symmetry are synonymous to 1/(1) x Z 2 
SPT phases of fermions. 

In the absence of interactions, these 2D TCIs are clas¬ 
sified by two integers Z 0 Z, the Chern number N and 
the mirror Chern number riM associated with occupied 
bands. Since the Chern number is defined without rely¬ 
ing on the mirror symmetry, for our purpose it suffices 
to consider systems with N = 0, for which the mirror 
Chern number um is defined as the Chern number of 
the occupied bands with the mirror eigenvalue +1 [32j . 
For example, (001) thin films of SnTe and monolayers of 
IV-VI semiconductors are predicted to be 2D TCIs with 
IumI = 2 [33H36]. 

To study the classification of 17(1) x Z 2 SPT phases 
in the presence of interactions, we follow the general ap¬ 
proach presented in the seminal work of Lu and Vish- 
wanath m and analyze the stability of noninteracting 
edge states against interactions. The existence of edge 
states that can only be gapped by breaking the mirror 
symmetry signals a 2D SPT phase. To begin with, the 
low-energy Hamiltonian for edge states of noninteracting 
TCIs is given by 

Ho = '^Vf / dx{-ill}\p.da:llla,R + ( 1 ) 

a 

Here the fermion fields '4’a,R/L denote respectively the a- 
th right and left movers (a = 1, ...,n), which transform 
differently under mirror: 


( 2 ) 
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where ry = sgn(nM)- The difference in mirror eigenvalues 
forbids single-particle backscattering between left and 
right movers; hence without interactions, gapless edge 
states are protected for any integer hm 0. The veloc¬ 
ity of different edge modes are chosen to be the same for 
simplicity; relaxing this condition will not affect any of 
our conclusions. 

We use bosonization to study the effect of interactions 
at the edge [3S1 [5H] • The bosonized Lagrangian for Hq 
takes the form 


L — . Kijdx4^idt^j . ; 

47r 47r 

(3) 

where K is an integer-valued matrix given by 


( Inxn 0 \ 

^^ 0 -1 

\ ^ -*-nxn/ 

(4) 

with Inxn being the n x n identity matrix, 
field (piix) satisfies the Kac-Moody algebra 

The boson 

[(j)i{x),da;'(j)j{x')] = 2mK~^5{x - x'] 

(5) 

and the fermion fields ' 0 ^ are given by 



(6) 


where different fields $are specified by a set of linearly 
independent integer-valued vectors La, a = To 

ensure that these fields can simultaneously have classical 
values, the commutator between any two of them must 
vanish. Since eq. <§ implies 

[^LSx),da,'^LA^')] =2TTiLlKLb5{x-x'), ( 11 ) 

this commutativity condition requires 

LlKLb = 0, ( 12 ) 


for any indices a, 5 = 1 ,..., n. A set of such vectors {La} 
will be referred to as a set of gapping vectors. As a 
general principle of bulk-boundary correspondence, the 
symmetry property of gapped edge states due to strong 
interactions reflects the topological property of the bulk. 
If the gapped edge preserves the U{l)xZ 2 symmetry, the 
bulk is in a trivial phase, i.e., adiabatically connected to 
an atomic insulator. 

We now show this scenario occurs for edge states that 
have n = 4 pairs of counter-propagating modes in the 
noninteracting limit. Such edge states can be gapped by 
interactions taking the bonsonized form eq. ( 10 ), with 
the following set of gapping vectors La- 


Electron interactions such as backscattering and 
umklapp processes can potentially gap the counter- 
propagating edge modes. These interaction terms are 
built from multi-electron creation and annihilation op¬ 
erators and are represented by cosine terms of the form 
cos{^l{x) + aL{x)), where the field ^l{x) = L'^K(j){x) 
is defined by an integer-valued vector L, and is an 
arbitrary phase. For our purpose, the interactions must 
preserve the charge conservation and mirror symmetry 
indispensable to 2D TCIs. It follows from eq. that 
charge conservation requires 

= 0, with t = (1„, 1„)^. (7) 

where 1 „ is the n-dimensional vector with all components 
equal to 1. For charge-conserving interactions, we further 
note the transformation law of the fermion field ([^ under 
mirror symmetry implies 

= ^l+ V^L’^rn, with m = ( 1 „, - 1 „)^. 

( 8 ) 

Hence the condition for mirror symmetry requires 

L^m = 0 mod 4. (9) 

To diagnose SPT phases, we consider sufficiently 
strong, symmetry-preserving interactions that com¬ 
pletely gap the 2n edge modes. This can be achieved by 
adding to the edge Lagrangian ^ n cosine terms [40] : 

n 

V = '^Xa C 0 S($I,„ {x)), (10) 

a=l 


Ti = ( 1 , 1 , 0 , 0 ;-l,-l, 0 , 0 )^, 

T2 = (0,0,1,1;0,0,-1,-!)^, 

T3 = (1,-1,0,0;0,0,-1,1)^, 

T 4 = (1,0,1,0;-l,0,-l,0)^. (13) 


It is easy to check that Li,...,L 4 satisfy the symmetry 
conditions 0 and (|^ , as well as the commutativity con¬ 
dition (HI). To motivate the choice of interactions ( [T^ , it 
is useful to regard four edge modes as two pairs of spin¬ 
ful Luttinger liquid in a two-leg fermion ladder system 
at half-filling. In the absence of inter-chain tunneling, 
the left- and right-moving modes have crystal momenta 
±7r/2 and transform oppositely under the lattice trans¬ 
lation: —>■ —ic\^. This is identical to the 

mirror symmetry transformation property of TCI edge 
states ([^— the only difference due to the factor i can 
be eliminated by redefining the symmetry operator [32] . 
Guided by this correspondence, we choose the interac¬ 
tions for n = 4 edge states denoted by Li and L 2 to be 
the bosonized form of the Hubbard interaction in the two- 
leg ladder, and L 3 and T 4 to be the antiferromagnetic 
inter-chain coupling. The former opens up charge gap 
and effectively generates two spin chains; the latter opens 
up a spin gap and leads to a rung-singlet phase that is 
fully gapped and translationally invariant. Equivalently, 
the interactions (13) gap the n = 4 edge states while 
preserving the mirror symmetry. A detailed analysis can 
be found in the Supplementary Material m- Therefore, 
we conclude that a noninteracting 2D TCI with mirror 
Chern number um = ±4 becomes trivial in the presence 
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of interactions. The additive nature of SPT phases then 
implies the same conclusion holds for um = 4A;, where k 
is an integer. 

Next we show case by case that the gapped edges of 
TCIs with n = 1 and 2 necessarily break the mirror sym¬ 
metry spontaneously. First, n = 1 edge states consist 
of a pair of counter-propagating modes, which can be 
gapped by symmetry-allowed Umklapp interactions that 
backscatter an even number of electrons from left to right 
movers, described by cos(2A:<i)£,) with L — (1,-1)^. The 
gap generation then implies 4)1, is pinned, i.e., ^ 0. 

This signals spontaneous mirror symmetry breaking, as 
can be seen from Q. 

For n = 2, by an exhaustive enumeration, we find 
two types of symmetry-preserving two-body interactions 
that gap the edge states, which are specified by two sets 
of gapping vectors {Li,L 2 \ and {Li,L 2 \ respectively, 
with Li = (1,1; —1, —1)^,il 2 = (Ij—1;—1)1)^ and 
L 2 = (1,—1;1,—1)^. We further note that the second 
type of interaction becomes equivalent to the first after 
redefining the flavor index of the left-movers 
Hence only the first type of interaction needs to be con¬ 
sidered. In terms of the electron operators, this interac¬ 
tion takes the form 

V =Ai(' 0 |^V'LV' 2 L' 0 il +h.c.) 

+ M('lplji'<plL'^lL'lp2R + h.C.). (14) 


Both terms conserve the number of fermions in each fla¬ 
vor (denoted by a = 1, 2) and commute with each other. 
The first term is an Umklapp process that backscatters 
two electrons with different flavors, and the second term 
flips the flavor of left and right movers simultaneously. 
It is convenient to introduce boson fields for each fla¬ 
vor: ipa = ((/>a,R + 4‘a,L)/‘2 and 6a = {(j)a,R “ (t)a,L)/‘^, 
with Ua = dxda being the density of electrons in flavor a. 


Equation (14) then becomes 


U = Ai cos(26»i -f 26*2) + A 2 cos(26»i - 202). (15) 


In the presence of this interaction, the edge becomes 
gapped when the fields 0i and 02 are both pinned. This 
leads to nonzero expectation values of single-particle 
backscattering operators: ~ 0 and 

(e* 2 ® 2 ) ^ (^/> 2 iiV’ 2 L) 7 ^ 0 , which implies spontaneous mir¬ 
ror symmetry breaking. 

The above edge state analysis shows that noninter¬ 
acting TCIs with mirror Chern number um = ±1 and 
±2 remain topologically nontrivial in the presence of in¬ 
teractions, contrary to the previous case of um = 4fc. 
Therefore, we conclude that interactions reduce the clas¬ 
sification of 2D TCIs protected by mirror symmetry, or 
U{1) X Z 2 SPT phases, from Z to Z 4 . 

In addition to its theoretical value, the above result has 
important implications for thin films/monolayers of SnTe 
and other IV-VI semiconductors, which are predicted to 


be 2D TCIs with {umI = 2 by band structure calcula- 
23H36] . Our analysis shows that interactions of the 


tions 

form (141 can qualitatively change the properties of n = 2 


edge states. At generic filling, only the flavor-flipping A 2 
term is allowed by momentum conservation and it is rel¬ 
evant for repulsive Luttinger interaction from the renor¬ 
malization group analysis m- As a result, there appears 
a gap in the flavor sector while the charge sector remains 
gapless and fluctuates. Boundaries and impurities affect 
the charge mode by pinning a fluctuating charge den¬ 
sity wave, which can be detected by STM measurement 
similar to the case of Luther-Emery liquid with a spin 
gap [42]. 

Interacting TCIs in three dimensions: We now 

turn to TCIs in three dimensions, protected by a single 
mirror symmetry, say x —>■ —x. Within band theory, one 
can define the mirror Chern number um on the 2D plane 
fca; = 0 in fc-space, which is invariant under this reflection. 
The integer um thus classifies 3D noninteracting TCIs [1] 
The hallmark surface states, present on crystal 
surfaces symmetric under mirror, consist of n = Itt-mI 
Dirac cones: 

n « 

Ho = ^vf / dr ipl{r){-idxSy + idySx)ilJa{r), (16) 


where is a two-component fermion field. 

Reflection acts on both electron’s coordinate and spin as 
follows: 


M'il>l{x,y)M ^ = Sxiljl{-x,y). (17) 

The mirror symmetry forbids any Dirac mass term 
and thus protects these n flavors of gapless Dirac 
fermions. 

Can the above Dirac fermion surface states be gapped 
by interactions without breaking the charge conservation 
and mirror symmetry? Finding the answer to this ques¬ 
tions will hold the key to the classification of interacting 
TCIs in three dimensions. This is a challenging task re¬ 
quiring non-perturbative approach to strongly interact¬ 
ing Dirac fermions in two dimensions. 

We now demonstrate explicitly that interactions can 
turn surface states with n = 8 flavors of Dirac fermions 
into a gapped and mirror symmetric phase without in¬ 
trinsic topological order (i.e., without fractional excita¬ 
tions). Such a completely trivial surface phase is con¬ 
structed as follows. First, we introduce a spatially alter¬ 
nating Dirac mass term to Hq: 



(18) 

where m{x) is a periodic function of x that alternates 
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FIG. 1: (Color online) Periodic array of ID domain wall 
fermions, generated by spatially alternating Dirac masses 
to 8 flavors of 2D Dirac fermions, see Eq.( |18|l^ . ID chi¬ 
ral fermion modes in flavors (red arrows) and flavors 

5,..., 8 (blue arrows) propagate in opposite directions along a 
domain wall. Counter-propagating chiral fermions have oppo¬ 
site mirror eigenvalues ±1. Each domain wall becomes gapped 
under the interaction ( fTopl ), thus leading to a gapped and 
mirror-symmetric 2D phase. 


(b) 


FIG. 2: (Color online) (a) A mass domain wall setup on 
a 3D TCI surface consisting of n flavors of Dirac fermions. 
One-dimensional chiral fermions reside at the domain wall at 
a; = 0, with n+ (n_ = n — n_) modes move in the +y {—y) 
direction, depending on the signs of Dirac masses mi,..., m„. 
(b) The domain wall in (a) is expanded to a wide region, 
sandwiched between semi-infinite regions on the left and on 
the right, with opposite Dirac masses. Importantly, (a) and 
(b) are both symmetric under mirror x —>■ —x and topologi¬ 
cally equivalent. For n 8k, the domain wall in (a), hence 
the middle region in (b) as well, cannot be gapped and mirror- 
symmetric. 



between toq and — toq, 


m{x) = 


mo for {2k — 1)L < x < 2kL, 
—mo for 2kL < x < (2k + 1)L. 


(19) 


Importantly, the resulting periodic array of Dirac 
mass domains preserves the mirror symmetry, because 
m{x) = —m{—x) and M'ij}\{x,y)sz'4’a{x,y)M~^ = 

-i’l{-x,y)s^il)a{-x,y). 

When the Dirac mass mo is large and the width of 
the domain L is large, the low-energy degrees of freedom 
are confined to the domain walls at a; = kL, where the 
Dirac mass changes sign. As is well-known, the mass do¬ 
main wall of a 2D Dirac fermion hosts a ID chiral fermion 
mode, whose directionality is reversed upon changing the 
signs of the Dirac masses on both sides. Therefore, our 
setup described by Ho + Hm hosts an array of ID domain 
wall fermions, one per flavor. On each domain wall, chi¬ 
ral fermions in flavors 1,..., 4 and those in flavors 5,..., 8 
move in opposite directions, and importantly, have op¬ 
posite mirror eigenvalues ±1 under the spatial reflection 
interchanging the two sides of the domain wall, as shown 
in Fig. [2 

We now draw a connection between the domain wall 
states on the surface of 3D TCIs to the edge states of 
2D TCIs: both are ID system of counter-propagating 
fermions with opposite mirror eigenvalues. Without in¬ 
teractions, the locking between the directionality and 
mirror eigenvalue forbids single-particle backscattering, 
leaving such ID system gapless. However, as we have 
shown earlier, the interaction given by eqs. (10) and (13) 


opens up a gap when there are four pairs of counter- 
propagating modes. Applying this interaction to each 
domain wall that we set up on the surface of noninter¬ 
acting TCIs then gaps the entire surface state states with 
n = 8 Dirac fermions, while preserving the mirror sym¬ 
metry X —)■ —x. We have thus explicitly constructed, us¬ 
ing a periodic array of domain walls, a completely trivial 


and gapped surface, the existence of which then implies 
noninteracting TCIs with mirror Chern number um = 8 k 
become trivial in the presence of interactions. 

Next, let us consider surface states of TCIs with um ^ 
8 k. Below we prove by contradiction that interactions 
cannot generate a gapped, mirror symmetric and non- 
fractionalized phase for these surface states [46]. Sup¬ 
pose such a trivial gapped phase exists, it must be adia- 
batically connectable to a massive Dirac fermion phase, 
where the Dirac masses are generated by external mirror 
symmetry breaking perturbations. This motivates us to 
consider a sandwich setup shown in Fig. [^, where this 
trivial phase takes up the region |a;| < L; to its right is a 
massive phase with a set of Dirac masses {wq}; and to 
its left the mirror image, a massive phase with opposite 
Dirac masses {—ma}. By construction, this sandwich 
setup is symmetric under the reflection x —)■ —x. 

We choose L to be much larger than the correlation 
length of the trivial gapped phase and let the surface 
Hamiltonian vary slowly with x across the interface at 
X = ±L, so that the trivial gapped phase (presumed to 
exit) adiabatically evolves into the massive Dirac fermion 
phase, without closing gap at the interface. Therefore, 
the surface is everywhere gapped and as a whole preserves 
the mirror symmetry. 

On the other hand, the sandwich setup is topologically 
equivalent to a domain wall between two domains with 
opposite Dirac masses (Fig. i^)- Without interactions, 
this domain wall hosts n = {umI flavors of ID chiral 
fermions, with n_|_ flavors and n_ flavors moving in op¬ 
posite directions and carrying opposite mirror eigenval¬ 
ues. Here n+ (n_) is the number of Dirac fermions with 
ma > 0 (toq < 0), and n+ + ri- = n. Importantly, for 
n+ ^ ri-, the domain wall must be gapless due to the 
presence of a net chirality, and for = n_ = n/2 ^ 4fc, 
we have shown earlier that the domain wall cannot be 
trivially gapped by interactions either. This result of the 
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domain wall contradicts with that of the sandwich setup, 
which is deduced to be gapped under the assumption that 
a trivial gapped surface is allowed on n ^ 8k TCI sur¬ 
faces. This contradiction proves the assumption wrong. 
Instead, 3D TCIs with mirror Chern number um ^ 8 k 
cannot have a trivial gapped surface and hence remain 
topologically nontrivial in the presence of interactions. 
Putting everything together, we conclude that interac¬ 
tions reduce the classification of 3D TCIs with mirror 
symmetry from Z to Z^. 

In addition to reducing the classification of nonin¬ 
teracting TCIs, interactions may also enable new TCI 
phases that do not exist in free fermion systems, as re¬ 
cently found in other symmetry classes [13148]. We leave 
this interesting problem of interaction-enabled TCIs with 
mirror symmetry for future study. 
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Supplemental Material 


In the Supplemental Material, we provide the approach 
to gapping vectors for n = 4 case, and the renormaliza¬ 
tion group (RG) analysis for n = 2 case. We utilize the 
Tomonaga-Luttinger liquid description in the following 
analyses, which we introduce first. 

In the Tomonaga-Luttinger theory, the fermion field 
'^l,R/L is given by [M] 


t _ ^ ^ ^i(‘Pa.+Sa) 

^I'Ka \/2TTa 

t _ ^ + ^ ^i{Va- 8 a) 

\j2'Ka \/2 'K(x 


( 20 ) 


a is an infinitesimal convergence factor. We neglect the 
Klein factor here. The two fields and daix) satisfy 

the commutation relation 


[0a{x), <^a{x)\ = i|sgn(a; - a;'), ( 21 ) 

and transform under the mirror symmetry as 

MiPa{x)M~^ = (fiaix), Mea{x)M~^ = 6a{x) + (22) 


The electron density for the right and left movers are 
given by ria^R/L = dx4>a(+n )/and thus the total elec¬ 
tron density of a-th pair is Ua = na,R + ria^L = dxOal'x. 
With the definition above, the bosonized Hamiltonian 
without gap-opening scatterings is 


H°=Y.Ha{Va,Ka) (23) 

a—1 


pairs of counter-propagating edge modes. It is instruc¬ 
tive to make an analogy between such edge states and 
the low-energy states of a two-leg fermion ladder system 
at half-filling. In the noninteracting limit, each chain is 
described by the tight-binding Hamiltoinan 

7Lo=-tE(SHi.-Gi<. + H.c.), (27) 

i.o- 

where j is a site index and the spin a =t, i- Each 
chain supports spin-degenerate left- and right-moving 
modes, which have crystal momenta -in: 12 respectively 
and transform oppositely under the lattice translation: 

(28) 

This is equivalent to the mirror symmetry transformation 
property of TCI edge states stated in the main text. 

We now add on-site Hubbard interaction to each chain, 

Hu = UJ 2 - 0 • (29) 

For [/ > 0, the repulsive interaction opens up a charge 
gap. As a result, at low energy each chain is equivalent to 
a spin-1/2 chain, as in the antiferromagnetic Heisenberg 
spin model. By further introducing antiferromagnetic 
inter-chain coupling, we obtain the rung-singlet phase of 
two coupled spin chains, which is gapped and transla- 
tionally invariant. Back to our original problem, this 
corresponds a gapped edge preserving the mirror sym¬ 
metry. 

The remaining task is to derive the bosonized form 
of the microscopic Hubbard and spin interactions in 
terms of the left- and right-moving fermion fields. The 
bosonized form of the Hubbard model is 


with 


2t: 


dx 


Kaidx^Paf + -^{dxOa 
Kn 


(24) 


The forward scattering terms, g 2 and < 74 , are included 
through the Luttinger parameter Ka and renormalized 
velocity Ua, defined as 


Va = VF,a 



1 + ( 94 ,a - 92,a)/{‘2T:VF,a) 

1 + (54.a + g2,a)/(2TrVF,a) ’ 


54 ,a 
2t:vf c 


f 92,a 
\2ttvfi 


(25) 

(26) 


Gapped states for n = 4 : rung-singlet phase 


H = Hp + H,, (30) 

where the “charge” (p) and “spin” (cr) degrees of freedom 
are separated to give 


2U r 

Hp = H°{vp,Kp) - — -^ / dxcos{2V20p), 

(27raj^ J 

2U f 

= H°(va,Ka) + — -^ / dxcos{2^/26a). 

(27ra)^ J 


(31) 

(32) 


The fields are defined by 

^p/o- = /^(^i ^ 2 ), V>p/cr = 7 ^ 2 ), (33) 


and the renormalized velocities and the Luttinger param¬ 
eters are 



1/2 

^ . U \ 

( ^ ^- 

II 

b 

1 ±- ) 

V 'XVF 2 

V 'XVF J 


- 1/2 


In this section, we elaborate on the interaction Hamil¬ 
tonian for n = 4 in the main text, which gaps four 


Xpjcr Xf 


( 34 ) 
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The Umklapp process in eq. (311 is relevant for Kp < 1 at 
half-filling, and the backscattering in eq. (32) is relevant 
for K„ < 1 at generic filling. For U > 0 and half-filling, 
the Umklapp process is relevant and opens the charge 
gap. This Umklapp interaction corresponds to two gap¬ 
ping vectors: 


Xi = (1,1,0, 0;-l,-1,0,0)^, 

X2 = (0,0,1,1; 0,0-1,-1)^. (35) 


Still the “spin” modes remain gapless. Now we can 
regard the two gapless modes as a two-leg ladder of spin 
chains [50] . The Jordan-Wigner transformation for a spin 
chain i (i = 1 , 2) and subsequent bosonization give 

S-{x) = -dj;0i(x) + ^ cos(26»j(j;)), 

TT Tra 

S+(x) = [(-!)- +cos(20,(x))]. (36) 

V ZiT a 

{x is actually defined on a lattice x = aj and finally 
we take the continuum limit. Here we take the lattice 
constant a = 1 for simplicity.) We assume the spin chains 
are written as the Heisenberg model 

= (37) 

3 

with i is the spin chain index and j denotes a site, i.e., a 
rung. The interchain coupling acts in a rung as 


H^=JiyY,{Sl^Sl^+Sl^Sy+JlY.Sl^SI^^. (38) 


The total Hamiltonian of the two-leg spin ladder is 

H = Hi + H 2 + (39) 

It can be decomposed as 

H^Hs + Ha, (40) 


where 



H, = H°ius,KY + 

to 

Ml 

f dxcos{2\/20s), 

(27rQ!)2 J 

Ha = H°iUa, Ka) + 

2 Jl 

j dx cos(2'/20a) 

(27ra)2 ^ 


+ 75—^ / f^a;cos(v2y)Q). (42) 

{ZTTaY J 

s denotes the symmetric part and a the antisymmetric 
part, defined by 

^s/a = + ^ 2 ), '^s/a ^ + l ^ + 2 )- (43) 

The renormalized velocity and Luttinger parameters are 


Here K = 1/2 since we assume the Heisenberg model 
for each spin chains. For the symmetric part (411, 
cos{2\/26s) is relevant for < 1 due to the RG analysis. 
In contrast, for the antisymmetric part ( |4^ , two cosine 
terms compete but from the RG analysis cos(2-\/20a) is 
relevant for Ka < 1/2 and cos(-\/2v^a) for Ka > 1/2. 
If we assume the antiferromagnetic interchain coupling, 
i.e., J/_ > 0, the Luttinger parameter becomes Ka> 1/2 
and thus cos(-\/2<7’o) is relevant. Now the two fields are 
pinned, the system becomes completely gapped. 

Since we assume that the interchain coupling is antifer¬ 
romagnetic, two spins in a rung form a singlet; it is called 
a rung-singlet phase. The way of gapping two-leg ladders 
is related to the Haldane gap for integer spin chains. 

The next step is to determine the corresponding gap¬ 
ping vectors. Note that the fields in the two-leg spin lad¬ 
der model come from the “spin” modes of the Hubbard 
model we considered first. Thus we should replace 


6*1 —>■ 6*1 — 6*2, 6*2 —t 03 ~ ^ 4 , 

^ — y>4, (45) 

to obtain L. By these replacement, we obtain the gap¬ 
ping vectors 

^3 = ( 1 ,- 1 , 1 ,- 1 ;- 1 , 1 ,- 1 , 1 )^, 

i:4 = (1,-1,-1,1; 1 ,- 1 ,-1,1)^, (46) 

where X 3 and L 4 correspond to cos( 2 -\/ 20 s) and 
cos(-\/ 2 (^a), respectively. 

Now we have four gapping vectors La (a = 1,...,4). 
However, we should confirm the absence of spontaneous 
symmetry breaking. An SPT phase and a trivial phase 
(such as an atomic insulator) are distinguished by the 
symmetry property of the gapped edge states. While a 
trivial phase permits a gapped and symmetry-preserving 
edge, edge states of a SPT phase, if gapped, must spon¬ 
taneously break the protecting symmetry. As shown by 
Levin and Stern ISD, spontaneous symmetry breaking 
may (but not necessarily) occur when a linear combina¬ 
tion of gapping vectors for the coefficients {ci} 

with no common divisors is nonprimitive, i.e., 

= cL (47) 

i 

and the integer c is larger than 1. In this case, the set 
of pinned fields {4>i.}, which themselves are symmetry¬ 
preserving, also freezes the field The latter may 

or may not break the original symmetry of the system, 
which needs to be checked case by case. Conversely, if 
CiLi is primitive for any coefficients with no common 
divisors, spontaneous symmetry breaking is guaranteed 
to be absent. 

The set of four gapping vectors Ta (a = 1, 4) is not 

primitive. Thus we define a new primitive set L/ from 













the linear combinations of La'. 

L' =i2 = (0,0,l,l;0,0,-l,-lf, 

L'g = \{L:, + Li) = (l,-l,0,0;0,0,-l,lf, 

L'^ = + X 2 + Lg) = (1,0,1,0; -1,0, -1,0)^. (48) 

The new set L'^ respects the 17(1) x Z 2 symmetry, and the 
primitivity ensures the absence of spontaneous symmetry 
breaking. Therefore we conclude that the edge modes can 
be gapped out without breaking symmetry for n = 4. 
The set of gapping vectors L'a is given in eq. (13) in the 
main text. 


RG analysis for n = 2 edge modes 

As mentioned in the main text, we find only two sets 
of symmetry-preserving gapping vectors for n = 2: 

Ti = (l,l;-1,-1)^, 

^2 = (1,-1;-1,1)^, (49) 

and 


The matrices and Mg can be diagonalized simulta¬ 
neously to obtain 

H J dx K+{da,^+f + 

+ ^ y da; K_{dx^p-f + f (56) 

with the new Luttinger parameter 



and the renormalized velocity 



The fields Lp± and 9± are defined by 

ip± = 9± = ^(9 i± 92). (59) 

First we consider the scattering processes denoted by 
Li and L 2 . The two scattering processes are written as 

Vi=gu J + li-c-), (60) 


Ti = (l,l;-1,-1)^, 

T2 = (1,-1; 1,-1)^. (50) 

In the following, we will consider the two sets in the 
Tomonaga-Luttinger description, and will analyze their 
relevance by RG calculation. At the beginning, we as¬ 
sume the two equivalent edge modes by setting vi = 
V 2 = V and Ki = K 2 = K. When two pairs of edge 
modes exist, two types of forward scatterings connecting 
two copies are allowed: 

12=52 J dx{i^lj^'lpiR-lpl^'lp2L +i’lL'^lL'lplR'^2R), (51) 
14=54 J dx{tplj^1piRlljlj^1p2R + 'tplL'^lL^lL'^2L). (52) 


Bosonizing the two processes 142 S'!!*! l 4 ^ obtain 

dx[{d,0fM^id,0) + {dJfMe{dJ)], (53) 


where 0 = {(pi,ip 2 )'^, 9 = (5i, ^ 2 )^, and the matrices 
and Mg are given by 


M^ = 
Mg = 


( vK (54 - 52)/27r\ 

U54 - 52)/27r vK ) ’ 


( v/K {g'^ + g^)/27r\ 
V(54 +52)/27r v/K )' 


(54) 

(55) 


I 2 = 56 J dx{\p\j^ilj\^il)iL'tl) 2 R + h.c.). (61) 

Vi is an Umklapp process occurring at half-filling kp = 
7r/2 and V 2 is a backscattering allowed at generic filling. 
Their bosonized forms are 

li = [ dxcos{29i + 292 ), (62) 

(27ra)^ J 

^2 = [ dxcos{29i - 292), (63) 

(27ra)2 J 

or by using 9± 

Vi = / dxcos{2V29+), (64) 

(27ra)2 J 

I 2 = / dxcos(2v^5_). (65) 

(27ra)2 J 

The RG analysis shows that Vi is relevant for 77+ < 1 
and V 2 for K- < 1. When a scattering process is rele¬ 
vant, it pins the field 9± and generates a gap. The pin¬ 
ning of 9± leads to the mass A+, estimated as A+ « 
(u+/Q;)(g„)^/^^“^*-+i and A_ « (u_/a)(( 7 h)^/i^“^*--i. 
This situation resembles the charge-spin separation of 
conventional spinful ID systems. The fields (pp and 5+ 
correspond to the charge degrees, and ip- and 5_ to the 
spin degrees. The charge sector is gapped by the Umk¬ 
lapp process and the spin sector by the backscattering 
process. 
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(a) 1L 2L 1R 2R (I’) 2L 1L 1R 2R 


FIG. 3 : Backscattering processes for n = 2 . There are 
two possible cosine terms that represent backscattering: (a) 
cos( 20 i —262) and (b) cos( 2 ^i—2^22). The energy dispersion of 
(a) might appear when two copies are related by time-reversal 
symmetry. When the velocity of two copies are different, the 
energy dispersion would be like (b). 

The mirror symmetry restricts the simultaneous gap 
opening of Ay and A_ because the pinning of 9 ± means 
the pinning of 9 i^ 2 - Since 9 i and 6*2 have a periodicity 
of TT, 9 i + 92 is pinned at either 0 or tt (mod 27 r) for 
Qu < 0 , and either 7r/2 or 37r/2 (mod 27 r) for > 0 . 
Similar consideration applies for the backscattering pro¬ 
cess, which pins 9 i — 6*2 and its value depends on the 
sign of 5b. Therefore 2 have expectation values of ei¬ 
ther 0 , 7r/4, 7r/2, or 37r/4, depending on the signs of gu 
and 5b, and hence the mirror symmetry is spontaneously 
broken. It is also leads to non-zero expectation values 
of single-particle backscattering ~ 7^ 0 

and ^ (^/'2 /jV'2l) 7^ 0, which is prohibited by the 

mirror symmetry. 

Next we consider Li and L2. L2 corresponds to 

V2 = gb J dx{^lj^'ijjl^i:2LilJ2R + h.c.), (66) 

and its bosonized form is 

^2 = [ dx cos(2(/?i - 2(^2) 

( 27 ra )2 J 

= (o ^\2 [ dxcos{ 2 V 2 (p_). 

( 27 ra )2 J 


V2 is equivalent to V2 by the redehnition > ip^L ^'^d 
'^2L '^iL ■ redefinition, the velocities of the two 

modes should be the same. When the two velocities are 
different, V2 and 1^2 read different scattering processes 

(Fig.i. 

Finally we extend the analysis to the case where the 
two velocities are different Vi ^ V2 as well as Ki ^ K2. In 
this case, the “charge” and “spin” degrees are no longer 
separated. Here we concentrate on ¥2- In the RG analy¬ 
sis, the scattering process V2 is relevant for 6 > 0 , where <5 
is a scaling dimension of 5b, i.e., the coupling constant gb 
transforms into X^gb under the scaling r = (x,t) —>■ Ar. 
The scaling dimension d is given by <5 = 2 -|- < 5 cos with 
Scos being a scaling dimension of cos( 2 (/ji — 2 <p 2 )- Fol¬ 
lowing Ref. [ 52 ], ^cos is calculated from the correlator 
K(r) = (cos[2(pi(r) - 2Lp2{r)] cos[2(^i(0) - 2(p2(0)]) as 

K{Xr) = A'^“=X(r). (68) 

If we assume an infinitely long system at zero tempera¬ 
ture, the Euclidean action after integrating 9 fields is 

where a 2 x 2 matrix L{q,u}) is defined as 

L{q,oj) = ( 70 ) 

TT 

Then the correlator K{r) will be 


K{r) = -e 


47(r) 


( 71 ) 


( 67 ) 


with 




/(r) = 


— f ( T -} 

27r / 27r' 


(Fii^ + ^22^ ^12 ^2l)- 


( 72 ) 


To perform the integrations over 5 and w, we differentiate J(r) with respect to x, and then impose the boundary 
condition /(O) = 0 . Following this procedure, we obtain 

B — Aq\ 


I{r) = 


451 ( 5 ! - 5 !) 

B — Ari 2 

452 (’ll - vf) 


log 


log 


a 

a — ix rjiT 
a 


log 


a — ix + 52T 


r 


log 


a 

a-\- ix -\- TfiT 
a 


a -\- ix -\- rj2T 


( 73 ) 


where 


A=^ + ^ + -{g'A + 92), 

Kl K2 TT 


B = (det Mg) 


viKi + V2K2 - -(54 - 92) 

TT 


( 74 ) 

( 75 ) 


’ll,2 = C T - (detM,p)(detMe), 


vl + vj 
^ 2 


(27r) 


2 ( 9 ?-9'i)- 


( 76 ) 

( 77 ) 























Note that 771^2 can be regarded as renormalized velocities. 
The scaling dimension Jcos becomes 
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Ariiri2 + B 
VlV2{'lll+V2)’ 


(78) 


and V 2 is relevant when S > 0, i.e., 

Ar]iV2 + B ^ 2 
1llV2{Vl+V2) 

For a simple case where 34 = 7/2 = 0; '^cos 


(79) 
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